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^.^ Abstract. The Alcubierre warp spacetime yields a fascinating chance for comfortable 

fi interstellar travel between arbitrary distant places without the time dilation effect as in special 

^H relativistic flights. Even though the warp spacetime needs exotic matter for its construction and 

I is thus far from being physically feasible, it offers a rich playground for studying geodesies in 

jH the general theory of relativity. This paper is addressed to graduate students who have finished 

a first course in general relativity to give them a deeper inside in the calculation of non-affinely 
parametrized null and time-like geodesies and a straightforward approach to determine the 
gravitational lensing effect due to curved spacetime by means of the Jacobi equation. Both 
topics are necessary for a thorough discussion of the visual effects as observed by a traveller 
inside the warp bubble or a person looking from outside. The visual effects of the traveller can 
be reproduced with an interactive Java application. 
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1. Introduction 

Even with the very latest propulsion technology we are far from being able to comfortably 
travel around in our own solar system, not to mention to travel to extra-solar planets. But 
even if one day we would have the necessary rocket drives and we could accelerate to nearly 
the speed of light to overcome the enormous distances, the time dilation effect of special 
relativity would be boon and bane together. In a certain way, general relativity could be 
a way out because it offers several mathematical solutions of the Einstein field equations 
like the Morris-Thorne wormhole [ 1 1 or the warp metric by Alcubierre [2 1 which could be 
used for interstellar travel. While a wormhole represents a shortcut in spacetime, the warp 
metric contracts and expands the spacetime locally and encompasses like a bubble a nearly 
flat region that is, in some sense, decoupled from the rest of the spacetime. Unfortunately, 
both mathematical solutions need exotic matter, which makes them physically unrealizable. 
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Nonetheless, these spacetimes offer a rich playground for studying geodesies in general 
relativity. While geodesies in the Morris-Thorne wormhole can be handled numerically in 
a quite straightforward manner, the time-dependent Alcubierre spacetime makes it necessary 
to transform the geodesic equation, the equation for the parallel transport of vectors along 
time-like geodesies, and the equations to determine the gravitational lensing effect caused by 
the warp bubble into their non-affinely parametrized form. Then, the resulting equations can 
be integrated numerically as usual. 

The aim of our article is to study in detail the influence of the warp bubble on null and 
time-like geodesies. For that, we first transform all the relevant equations mentioned above 
into their non-affinely parametrized form. Beside the paths of light rays, we also discuss the 
frequency shift and the effect on a bundle of light rays resulting in gravitational lensing that 
let us visualize point-like objects in a most realistic way. By means of the parallel-transport 
equation, we show how particles initially at rest will be carried along and undergo geodesic 
precession. Furthermore, we explain what an observer would actually see inside or outside 
the warp bubble using either four-dimensional ray tracing or an interactive Java application. 
For the latter, we make use of the highly efficient hardware architecture of programmable 
graphical processing units (GPUs) to correctly visualize point-like stars in contrast to the 
elongated stars shown in science fiction movies. The view from inside the warp bubble will 
be compared to the view of a special relativistic observer. The differences between the views 
of these observers can be clearly reproduced using our Java application. 

A detailed discussion of null geodesies from inside the warp bubble, in particular the 
emergence of an apparent horizon behind the warp bubble, was given by Clark et al. [3 1 What 
an observer within the Alcubierre spacetime would actually see was shown by Weiskopf |4|. 
There are several papers that discuss the physical nature of warp drive spacetimes, and we can 
give only a few references. 0|6l|7]|8) 

To study geodesies in detail, we refer the reader to the interactive visualization tool 
GeodesicViewer ll9l [T0ll and how it could be used in the classroom. 

The structure of this paper is as follows. In section |2j we briefly review the Alcubierre 
warp metric and present the local reference frames of a comoving and a static observer. 
In section [3] we study the trajectories of light rays and how they influence the view of an 
observer that either co-moves with the warp bubble or stays static in the outside. We also 
discuss the frequency shift and the lensing effect caused by the warp bubble. For the case of 
a comoving observer, we present our interactive Java application and the necessary technical 
and implementation details in section HI Finally, we discuss the influence of the warp bubble 
on particle trajectories in sectionB] The technical details of the integration of the non-affinely 
parametrized equations are relegated to the appendix. 

The Java application and its sources, high-resolution images of this article, as well as 
some movies that show the visual distortion of the warp bubble or the motion of particles can 
be downloaded from http://www.vis.uni-stuttgart.de/~muelleta/Warp 

2. Warp metric and local reference frames 

The Warp metric developed by Miguel Alcubierre [2] can be described by the line element 

ds 2 = -c 2 dt 2 + [dx-vf(r)dt} 2 +dy 2 + dz 2 , (1) 



where c is the speed of light, 

dx(t) 
~ dt '' 



(2) 
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r(t) = ^(x-x(t)y+yl+z\ (3) 

_ tanh(g(r + fl))-tanh(g(r-fl)) 
J{ ' 2tanh(atf) ' ( ' 

and jc(f ) is the worldline of the center of the warp bubble. The parameters R > and a > in 
the shape function f{r) define the radius and the thickness of the bubble, see figurefT] 




Figure 1. Shape function f(r) for a radius R = 2 and thickness parameters o~ = 1 (solid line) 
and o~ = 8 (dashed line). 

For this metric, we can define two natural local tetrads en\ — <?^<9 U that represent the 
local reference frames of either a comoving or a static observer. The comoving tetrad is 
defined by 



■5(0) 

and the static local tetrad reads 



e (0) = ~(dt + vfd x ) , em = d x , e (2 ) = d y , e (3) = d z , 



1 



e (0) 



J* 



- - v 2 f 2 

yf 



d t , 



'(2) 



dy, 



-(3) 



c\[c<- 



-A 



Vc 2 ~v 2 f 2 



dx. 



(5) 



(6) 



(7) 



A _ v 2p 

It is obvious that the comoving tetrad is valid everywhere, whereas the static tetrad is defined 
only in the region of the spacetime where v 2 / 2 < c 2 . Both tetrads fulfill the orthonormality 
condition g ^e^ej ., — I7(;)(/) with T\t{\(f\ = diag(— 1, 1, 1, 1), which means that these tetrads 
locally define a Minkowskian system. 

Throughout the paper we consider a warp bubble that moves with constant velocity. 
Thus, the center of the warp bubble follows the worldline x(t) = vt. Furthermore, we set 
{c} = 1 for numerical examples. Then, times and distances are given in seconds and light- 
seconds or years and light-years, respectively. 



3. Null geodesies 

In this section, we will discuss the influence of the warp bubble on the propagation of light 
for several different situations. In general, light paths follow from the numerical integration 
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of the null geodesic equation 

£/V n dx V d X P _ 

~d?J +rv Plx^x-°' (8) 

where X is an affine parameter and Fy P are the Christoffel symbols of the second kind (see 



App. Appendix A for the Christoffel symbols of the Warp metric). After each integration step, 



we have to check that the constraint equation 
dx* dx v 
8 ^lXlX ={) (9) 

with metric tensor g^ v is still fulfilled. This constraint ensures that the geodesic remains 
light-like even if numerical error is present. 

For the warp metric, an integrator with step-size control is indispensable. Otherwise, the 
step-size would have to be inefficiently tiny. However, irrespective of the numerical integrator, 
direct integration of the geodesic equation leads to numerical problems at the rim of the bubble 
for certain initial values because of the inappropriate affine parameter. Although the spacetime 
coordinates (t,x,y,z) are smooth, the step-size of the affine parameter becomes extremely 
small. Thus, an integrator with step-size control will get stuck and the constraint equation will 
be violated. To avoid the numerical difficulties, we use the non-affinely parametrized geodesic 



as described in App. Appendix B Additionally, the resulting equations are numerically much 
more accurate than the ones that follow from the affinely parametrized geodesic equation. 

The gravitational frequency shift Zf between the emitted, (Qsrc, and the observed, O ^ s , 
light frequencies is obtained by 

, ®src gnvU S rc K , im 

w obs ^v"o bs ^ 
where Usrc is the four- velocity of the light source and u ^ s is the four- velocity of the observer, 
see e.g. Wald [ 1 1 1 . The tangent of the light ray is given by W" = dx^ jdX and must be evaluated 
either at the source or the observer position. If — 1 < if < 0, we call it a blueshift, and if Zf > 0, 
it is a redshift. 

To determine the lensing effect caused by the warp bubble, we study the behavior of 
the spacetime curvature on a bundle of light rays that is described by two Jacobian fields 
Y, = Yf <9jj . The change of the Jacobian fields along the central light ray with tangent k = W 1 d^ 
is determined by the Jacobian equation 

D 2 T M 

-^=R» vpa k v *PYf. (11) 

The cross section of the light bundle follows from the projection of the Jacobian fields onto 
the parallel-transported Sachs vectors s, = sf d^ that are perpendicular to the light ray k. The 
resulting Jacobi matrix 

Jtl = gwY?s) (12) 

describes how the shape of the initially circular bundle of light rays transforms into an ellipse 
with major and minor axes a± along the central light ray. Fortunately, we can calculate the 
matrix by first integrating from the observer to the source and then inverting this matrix to 
obtain the behavior of a light bundle from the source to the observer. The result of this 
calculation is the magnification factor /Zmag = A 2 /(a + a_). Details of this calculation in the 



context of the non-affinely parametrized geodesic equation can be found in App. Appendix 
[P] A thorough discussion of gravitational lensing, however, is out of the scope of this article 
and we refer the interested reader to the standard literature [ 12 1 (See also Ohanian fTJI for a 
discussion of the gravitational lensing by a Schwarzschild black hole.) 
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3.1. View from the bridge 

The view from the bridge was already discussed in detail by Clark et al. Here, we will 
reconstruct their results for the sake of completeness. Additionally, we will show how light 
rays behave in the close neighborhood of the warp bubble and how the warp bubble influences 
the lensing of point-like objects. 

The local reference frame of the bridge is given by the comoving tetrad of Eq. Q, 
Because of the axial symmetry, we can restrict to geodesies in the xy-plane. Then, an incident 
light ray with angle ^ with respect to the local reference frame can be described by the four 
vector 



0)( 



-cos^e (1) +sin^e( 2 )) 



= *%. 



(13) 



-e(o)-t-cosL,e (1 )- 

Here, we use the minus sign in front of the time-like tetrad vector em) because we integrate 
geodesies back in time. If we are only interested in the paths of the light rays, we can set the 
frequency co = 1 . 

Figure|2]shows several light rays for an observer on the bridge when the observer passes 
the origin at t = 0. Although all the incident light rays are equidistantly separated by 10° with 
respect to the observer's local reference frame, they apparently approach the observer from 
behind. 




Figure 2. Light rays reach the comoving observer on the bridge, here at (t = 0,.v = 0,y = 0), 
with incident angles if; = {10°, 20°, . . . , 170°} with respect to the comoving local reference 
frame eyi . The parameters of the warp metric read R = 2, o~ = 1, v — c. 



Figure [3] shows the view from the bridge of an observer moving with different warp 
speeds when passing the origin x = y = at t = 0. To visualize distortion effects close to the 
warp bubble, we use two checkered balls of radius r^ji = 0.5. The green ball is located at 
x = 10, y = whereas the red ball is located at x = 0,y = 3. Similar to the special relativistic 
motion in flat Minkowskian spacetime, there is an aberration in the direction of motion that 
is stronger the faster the warp bubble moves. In contrast to the special relativistic motion, 
however, the aberration here is only due to the curved spacetime. 

As already found by Clark et al. there is no aberration for £, = 90°. Hence, light rays 
that originate precisely at 90° to the direction of motion of the warp bubble will be seen 
by an observer at the bridge at an angle £, = 90° (see also App. Appendix C I. Furthermore, 
the observer on the bridge cannot see the whole spacetime. If we trace light rays from the 
bridge with initial angles t, back in time until they hit a sphere in the asymptotic background, 
t = /"max, cf. figure HI we find that there is an apparent horizon opposite to the direction of 
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Figure 3. View from the bridge at t = in the direction of motion for velocities v = 0.01c (a), 
V = c (b), v = 2c (c) and a panorama camera with 180° x 60° field of view. The Milky Way 
background is represented by a sphere with radius r ma x = 200. fiTi 



motion, where the size of the apparent horizon depends on the velocity of the warp bubble. 
However, the observer will not see a black region because for any direction % some region of 
the asymptotic background will be visible. 




20 40 60 80 i 100 120 140 160 18 



Figure 4. Angle <p at r ma x = 5 x 10 4 with respect to the incident angle £ for velocities 
v/c = {0.5, 1,2,9} and warp parameters R = 2, o~ = 1. The angle | = 0° corresponds to the 
direction of motion. 



To determine the frequency shift and the lensing caused by the warp bubble, we use the 
asymptotic sphere as the place where there are everywhere point-like light sources with unit 
frequency «src = 1- Then, the observer at the bridge will detect frequency shifts as shown in 
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figure p] Irrespective of the velocity, there is no frequency shift for t, = 90°. In the direction 
of motion, | < 90°, we have a slight blueshift, whereas in the opposite direction the redshift 
is quite strong. 
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Figure 5. Frequency shift 1 + Zf between r max = 5 x 10 4 and the bridge with respect to the 
incident angle t, for velocities v/c = {0.5, 1,2,9} and warp parameters R = 2, a = 1. 

The lensing of point-like sources caused by the warp bubble is shown in figure [6] In 
the direction of motion, we have a slight magnification, whereas in the opposite direction, the 
light of point-like sources are strongly dimmed. 
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Figure 6. Lensing /imag between r max = 5 x 10 4 and the bridge with respect to the incident 
angle t, for velocities v/c = {0.5, 1,2,9} and warp parameters R = 2, a = 1. 



Hence, even though the horizon does not produce a black region, the strong redshift 
together with the magnification let the region | < 90° appear dark for high velocities. 

The 'view from the bridge' for a special-relativistically moving observer is well- 
known from standard literature. The aberration and Doppler-shift can be derived from the 
representation of an initial direction k with respect to either the moving reference frame 

5(0) = ~(d t + v SI d x ) , e (1) = 7 ( v sr - + d x ) , (14) 
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= d y , 5(3) = d z , as in Eq. (13 1, or the standard Minkowskian frame. Here, vsr is the 



e( 2 ) = o y , e( 3 ) 

special-relativistic velocity of the moving observer. Thus, we obtain the aberration formula 

1 
D 



cos<r = ^(cos|-v sr / c ) (15) 



and the Doppler-shift 

l+z f = D, (16) 



where D = y[\ — (vsr/c) cos ^] is the Doppler-factor and y= 1/y 1 — v\ r /c 2 . Here, the primed 
angle is with respect to the flat Minkowski spacetime and the barred angle is with respect to 
the special-relativistic reference frame. The magnification factor is given by 

Mmag=£>~ 2 , (17) 

see Weiskopf et al. [ 14 1 for a detailed discussion. For velocities vsr close to the speed of light, 
we have very strong blueshift and huge magnification in the direction of motion. 

While the direction of no frequency shift, Zf — 0, in the Warp metric is fixed by 
E, = 90°, this borderline is defined by D = 1 in the special-relativistic case, which also 
gives the borderline of no-magnification, /imag = 1- The corresponding angle reads £, = 
arccos [(7— l)/(y/3)], where £, < 90° for vsr > 0. Hence, for velocities close to the speed of 
light, a special-relativistic traveller will only see a small bright spot in contrast to a traveller 
inside a warp bubble. 

3.2. Warp bubble approaching a static observer 

Consider a static observer located at x = 10, y = and a warp bubble with R — 2 and a = 1 that 
approaches the static observer with v = 2c. At t = 0, the bubble crosses the origin x = 0,y = 0. 
The observer, however, that looks towards the approaching warp bubble, will not recognize 
any significant image distortion, see figure [71 a). At t = 3.5 (figure I7tb)), the warp bubble is 
located at x^ — 7 and, thus, is already quite close to the observer that will see that some portion 
of the sky apparently shrinks in the direction of view. However, in the next few instances the 
view of the observer changes dramatically. Some parts of the sky in the direction of view 
disappear as a result of an apparent horizon like in the bridge observer example, see also 
Figs.|8]and|9] 



The corresponding frequency shift and lensing diagrams are shown in Figs. 10 and 11 
Here, ^ = 180° represents the direction to the approaching warp bubble. For t — 0, there is no 
significant frequency shift or magnification in any direction. For t = 3.98 the warp bubble is 
located at xu = 7.96, however, the approaching warp bubble leads to a strong blueshift in the 
direction of the warp bubble and a redshift in the opposite direction. The lensing effect lets 
stars appear brighter especially in the direction of the approaching warp bubble. 

3.3. Warp bubble passing the observer 

Consider a static observer located at x = 0,y = —4 and a warp bubble with parameters R = 2, 
a = 1 that passes the observer with v = 2c. At t = 0, the bubble crosses the origin x = 0,y = 0. 



Figures 12 and 13 show several representing light rays that reach the static observer at times 
t = 3 and f = 6 with incident angles ^ = {0°, 10°, ... , 180°}. The corresponding views for a 
panorama camera with 180° x 60° field of view heading towards the origin are shown in figure 



14 Here, the Milky Way background is represented by a sphere with radius r m ax = 200, and 



the checkered ball of radius r^ a jj = 0.5 is located at x — Q,y — 3. 
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Figure 7. View of a static observer located at x — 10, y = in the negative x-direction 
towards the approaching warp bubble that moves with velocity v = 2c and that crosses 
the origin x = y = at t = 0. The panorama camera has 180° x 60° field of view. The 
observation times are t = {0,3.5,3.98} when the bubble is at.*;, = {0,7,7.96}. The red ball 
is located at x = 0,y = 3. The Milky Way background is represented by a sphere with radius 
rmax = 200.E2] 




Figure 8. Incoming light rays for a static observer at (t ■ 
parameters R = 2, o~ = 1 , and v = 2c. 



10, y = 0) and warp 



The three phantom images of the ball (figure fl4[c)) follow from light rays with 
incident angles \ sa {28. 8°, 45°, 83.2°}. The distortion of the Milky Way background can 
be understood by means of the relation between the intersection angles (p of the light ray 
with the background sphere and the incident angles £,, see figure 15 The discontinuities 
reflect the apparent horizons at <i; nor s 
observation times t 



{171.3°, 102.8°, 153.6°, 164.8°} for the corresponding 
{0,3,6,9}. Since the mapping <p t-> § is non-injective, some parts of 
the sky appears more than once. For example, the point (r = r m axj<p = 50°) appears three 
times at observation time t = 9 under the incident angles t, w {10.82°, 30. 83°, 67. 85°}. 
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Figure 9. Angle (p at ("max = 5 x 10 4 with respect to the incident angle | for a static observer 
at x = 10, y = and t = {0, 3.5, 3.98}. The warp parameters read R = 2, (7 = 1, and v = 2c. 
The angle £ = 180° corresponds to the direction of the approaching warp bubble. 
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Figure 10. Frequency shift 1 + Zf between r ma x = 5 x 1 4 and a static observer at x = 1 0, y = 
and t = {0, 3.5, 3.98}. The warp parameters read R = 2, a = 1, and v = 2c. 



What can also be read from figure 15 is that, after the warp bubble has passed the 
observer, there is another distortion region that apparently moves in the negative ^-direction. 
This secondary distortion region can be easily understood. Due to the finite speed of light, 
light rays that have already traversed the warp bubble region at earlier times now reach the 
observer. Thus, by increasing time, the observer receives ever earlier light rays. 

Figure 16 shows the frequency shift for the observation times t = {0,3,6,9}. At t = 0, 



the frequency shift is negligible. However, when the warp bubble has passed the origin, 
there is some blue- and redshift from the primary distortion region (t, > 90°), whereas in the 
secondary region there is a strong redshift close to the apparent horizon directions. 



Figure 17 shows the lensing effect for the observation times t = {0,3,6,9}. Similar 
to the frequency shift, the magnification at t = is negligible. However, when the warp 
bubble passes the observer, there are strong magnifications at both sides of each bubble rim. 
In between, stars would appear dimmed. In the regions where the redshift is dominant, the 
lensing effects fades out the star light. 
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Figure 11. Lensing £lmag between r max = 5 x 10 4 and a static observer at x = 10, y = and 
t = {0,3.5,3.98}. The warp parameters read R = 2, a = 1, and v = 1c. 




Figure 12. Incoming light rays for a static observer at (t = 3,x = 0,y = —4) with incident 
angles b, = {0°, 10°, . . . , 180°}. The warp parameters read R = 2, CT = 1, and v = 1c. The 
colored lines are only for better distinguishability. 



4. Interactive flight 

In contrast to the usual graphical representation of the visual appearance of a star field in 
science fiction movies with its elongated stars, the actual view of such a star field as seen 
from within a warp bubble is completely different. To demonstrate what could be really seen, 
we have developed a Java application with which the user can travel through the stars of 
the Hipparcos catalogue either with relativistic speeds in the Minkowskian spacetime or with 
warp speed in the Alcubierre metric. 

In the following subsection, we briefly describe the relativistic star flight simulator 
JRelStarFlight and its graphical user interface. Technical details about the implementation 
and the calculation of the physical quantities are discussed in subsection|4.2| 
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Figure 13. Incoming light rays for a static observer at (t = d,x = Q.y = —4) with incident 
angles S, = {0°, 10°, . . . , 180°}. The warp parameters read R = 2, a = 1, and v = 2c. The 
colored lines are only for better distinguishability. 



4.1. JRelStarFlight 

The graphical user interface of our relativistic flight simulator JRelStarFlight is shown 
in figure [18] When the application is launched, the observer is at rest in a Minkowskian 
spacetime and uses a full-sky panorama camera. To toggle between special relativistic and 
warp flight, use the drop-down menu in the lower left corner. There are two camera modes: 
the full-sky camera maps the whole sky from spherical coordinates onto a rectilinear like grid; 
the pinhole camera acts as a normal camera with 50° vertical field of view. 

Moving the mouse in the OpenGL window with the left mouse button pressed, the 
viewing direction and, thus, the direction of motion can be changed. The current velocity 
j3 = v/c can be modified by means of the input field 'beta'. Toggling the 'play' button lets 
the observer move or stop. The current position can also be set with the input field 'distance 
to origin' . Please note that when the observer has moved away from the origin, the mouse 
rotation acts as if the observer has moved in the newly selected direction ab initio. 

4.2. Technical details 

The apparent position of a point-like star as seen within the warp bubble depends on how a 
light ray of this star is influenced by the curved spacetime. For the visual appearance, we also 
have to take the frequency shift and the lensing effect into account. 

For the actual visualization of the stars we adopt the method by Miiller and 
Weiskopf [18]. As database we use the Hipparcos star catalogue [19| and assign a Planck 
spectrum to each star with temperature r s t ar calculated by Reed's empirical law [20 1. The 
frequency-shifted spectrum is again a Planck spectrum but at the different temperature 
^obs = ^star/(l + £/)■ The resulting apparent visual magnitude of a star is determined by the 
magnification factor /imag and the integral over the Planck spectrum in the visual wavelength 
domain. A finite eye or telescope aperture yields a Fraunhofer diffraction pattern. 

Since the star distances are in terms of light-years, we use years as unit of time. With 
25 frames per second and a step size of At = QAy per frame, one second of visualization time 
corresponds to 10 years of simulation time. 
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Figure 14. View of a static observer located at x = 0,y = —4 in the positive y-direction for 
t = (a), f = 3 (b), t = 6 (c), and t = 9 (d). The warp parameters read R = 2, o~ = 1, and 
V = 2c. The panorama camera has 180° x 60° field of view. The ball is located at x = 0,v = 3. 
The Milky Way background is represented by a sphere with radius r^ = 200. 



4.3. Implementation details 

Before running the Java application we have to generate the two-dimensional lookup table 
that stores the observation angle ^, the frequency shift z/, and the magnification factor /Xmag 
for each asymptotic light direction <p and velocity v of the warp bubble. For that, we integrate 
the geodesic equation, the parallel transport of the Sachs basis vectors, and the Jacobian 
equations from the observer back in time until the geodesic reaches the sphere with radius 
'"max = 5 x 10 4 . Hence, we obtain a relation between observation angles | and intersection 
angles <p with the sphere, which must be inverted at the end of the calculations. 

For our simulator, we use the warp parameters R = 2, a = 1 and a lookup table with 
resolution 4096 x 248. While the angle <p is sampled linearly from to %, the velocity in the 
row number 'row' is given by 

v=10 (row+i)/248_ 1 (18) 



Thus, the minimum velocity v mm 
v'max = 9 in row = 247. 



! 0.0093 is stored in row = and the maximum velocity 
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Figure 15. Angle (p at r ma x = 200 with respect to incident angle t, for / = {0,3,6, 9} and an 
observer located at x = 0, v = —4. The warp parameters read R = 2, (7 = 1, and v = 2c. 
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Figure 16. Frequency shift 1 + z,f between r max = 200 and the observer located at x = 0,_v = 
—4 with respect to incident angle £, fort = {0,3,6,9}. The warp parameters read R = 2, (7 = 1, 
and v = 2c. 



The rendering of the stars is realized using a Java implementation of the open graphics 
library OpenGL and the shading language GLSL [21 1. In the vertex shader, we use the 
above explained lookup table to calculate the apparent positions and frequency-shifted Planck 
temperatures of the stars. In the fragment shader, the Planck temperatures are mapped to the 
precalculated color and Fraunhofer diffraction pattern. 



5. Timelike geodesies 

So far, we considered only the influence of the warp metric on light rays. Equally interesting, 
however, is the influence on massive particles as we will discuss in this section. As with the 
null geodesies, we use the non-affinely parametrized geodesic equation with the coordinate 
time as parameter. 
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Figure 17. Magnification ,Umag between r ma x = 200 and the observer located at.v = 0,y = —4 
with respect to incident angle \ for t = {0,3, 6, 9}. The warp parameters read R = 2, a = 1 , 
and v = 2c. 
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Figure 18. Screenshot of the Java application JRelStarFHght. 



5.1. Particles from the bridge 

A particle at the center of the warp bubble with zero initial velocity will stay there for ever 
irrespective of the worldline x(t) of the bubble. So, let us consider particles with initial 
local velocity Vp art with respect to the comoving reference frame and the corresponding four- 
velocity 

u = cy [e (0) +J5 (cos£e (1) +sin^e (2) )] = w%, 



(19) 



19 



shows particles that were emitted radially 
= 0.5c. Each solid line represents all 



where j3 = Vp all /c and y = l/-y/l — J3 2 . Figure 

from the bridge with initial local velocity Vp^ 

particles that have travelled a specific coordinate time Af . The dashed lines represent particle 

trajectories for a few initial directions ^ . 

In the first few seconds, the particles still move within the warp bubble and are carried 
along with it. Depending on their initial direction ^, they can leave the warp bubble after 
some time At. 

To measure the current velocity of a particle during this period, we can use only a valid 
observer, which is an observer represented by a comoving local tetrad. Such an observer will 
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Figure 19. Particles emitted radially from the bridge with initial local velocity i'p ar [ = 0.5c 
after Af = {1.0, 1.5, . . . , 6.0} (solid curves). The dashed lines correspond to particle trajectories 
for initial directions t, = {0°,30°,. ... 150°}. The warp parameters readi? = 2, £7 = 1, v = 2c. 



always measure a velocity less than the speed of light. Figure 20 shows the velocities of 
the particles after Af = 1000 when all of them are far away from the sphere of influence of 
the warp bubble. Similar to null geodesies, a particle with initial direction ^ = 90° is only 
displaced by the warp bubble and the velocity keeps unchanged. All other particles have either 
a higher or a lower velocity depending on the initial angle ^ . 
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Figure 20. Particle velocities Vp ar [(r = 1000)/c far outside the bubble's sphere of influence 
after Af that were emitted radially from the bridge with initial local velocity Vp ar t(r = 0) = 
0.5c. Particles with initial angles b, < 5° are still in the sphere of influence of the warp bubble. 



5.2. Particles injected from outside 

Consider a static observer located at x — Q,y — —4 that emits particles with local velocity 
v part = 0-5c in the directions ^ <G [0°, 180°] at coordinate time t. The corresponding four- 
velocity u with respect to the static reference frame reads 

-j3(cos§e (1) + sm§8( 2 ))]=w' 1 5 M) (20) 

The solid lines of figures 2 1 and 



u = cy [e (0) - 

where j3 = vp^fc and y = l/^/l — j8 
that were emitted at coordinate times t — — 6 or t = 
some time At. The thick solid lines correspond to t 
trajectories for a few initial directions £, . 



t-l 



22 



show particles 
4, respectively, and that have travelled 
0. The dashed lines represent particle 
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Depending on the emission time, the particles reach the sphere of influence of the warp 
bubble under different impact angles. Hence, they will be deflected, carried along, and 
accelerated or decelerated in different ways. 




Figure 21. The solid lines represent particles that start with fi = 0.5 from the position 
(t = —6,x = 0,y = —4). The lines are separated by A/ : = 0.5. The thick solid line indicates 
t = 0. The warp parameters read R = 2, a = 1, v = 2c. The dashed lines correspond to timelike 
geodesies with initial angles t, — 0° , 30° , . . . , 1 80° . 




Figure 22. The same situation as in figure|2ljbut with initial position it = — A,x = 0,y = —4). 



5.3. Particle field initially at rest 

Consider a field of 5 x 6 = 30 particles represented by tiny balls that are initially at rest with 
axes in the positive x direction. The initial position of particle (i,j) is given by x, = — 3 + i, 
yj — —4 + j. At the beginning of the simulation, t = —3, the warp bubble with parameters 
R = 2, a = l,v = 2c is located at x = — 6,y = 0. To determine the influence of the warp 
bubble on this particle field, we calculate the parallel transport of each particle. For that, 
we have to integrate the geodesic equation with initial conditions x^\ „ = (— 3,Jc,,y/,0), 
dx^ /dx\ „ = cem) together with the parallel-transport equation 

dp" M dx v dp^_ 

dx + vp dx dx ' 
where the four-vector p^\ „ = (0,1,0,0) describes the orientation of the particle. Like 
the geodesic equation, we first convert the parallel-transport equation into its non-affinely 
parametrized form (compare the parallel transport of the Sachs basis vectors explained in 



(21) 



App. Appendix D i 
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When the warp bubble approaches the particle field, the particles will be carried along 



for some time depending on the distance to the center of the warp bubble, see figure 23 

Additionally, the particles undergo a geodesic precession away from the jc-axis while they 
are in the sphere of influence of the warp bubble. The precession angle is larger the closer the 
particle is to the rim of the warp bubble. Particles at y — z — keep their orientation. 

6. Summary 

The Alcubierre warp spacetime offers a rich playground for studying geodesies in general 
relativity. For a direct numerical integration of the geodesic equation, the parallel transport, 
and the Jacobian equation, we have transformed these equations into their non-affinely 
parametrized form. 

The view of an observer comoving with the warp bubble is distorted similarly to the 
special relativistic aberration. In the direction of motion, however, the distortion is much 
softer than the special relativistic aberration. Additionally, for velocities larger than the speed 
of light, there is an event horizon behind the warp bubble. However, there is no black region 
as in the Schwarzschild spacetime. Nevertheless, the strong redshift as well as the large 
attenuation let the rear side appear very dark. 

The view of a static observer outside the warp bubble heavily depends on his location 
and observation time. Light rays are strongly deflected and multiple images might appear. 

An interactive special-relativistic flight or warp-drive voyage through the stars is 
simulated within our Java application. In contrast to the elongated stars shown in science 
fiction movies, stars would keep their point-like shape. 

Particles that start from inside the warp bubble, will be 'decelerated' or 'accelerated' 
depending on their initial angle to the warp bubble's direction of motion. Additionally, 
particles undergo geodesic precession while they are in the sphere of influence of the warp 
bubble. 

Appendix A. Christoffel symbols 

The Christoffel symbols of the Alcubierre metric (fill read 

K = f -^-, r-;; = -// y v 2 , r~, = -//y, r* = ^4^-/Av 2 -/ r v-/«9,v, (ad 



r t /fa 2 r x _ ffxV* r y _& r z _ & rA ™ 

c A c A 2 2 

r , _ _ffyf_ r * _ / 2 / v v 3 +c 2 / y v _ _.ff^ _ f 2 f z v 3 +c 2 f z v 

2c 2 ' l ty~ 2c 2 » l tz 2c 2 - ift" 2c 2 . V^*) 

pt _ fxV x _ ff.y -t _ ty x ffyV 2 

xx c 2 i L xx c 2 ' L xy 2( .2 ' L xy 2( p. ' V^-^J 

p _& r _ f/y (A 5) 

xz 2c 2 ' z 2c 2 ' v^- J ' 

with derivatives 

f = d f( r ) = -Hx-x s (t))df{r) 
dt r dr 

_df{r) _ x-x s (t)df{r) 
dx r dr 
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Figure 23. Parallel transport of tiny balls that are initially (t = —3) at rest. The dashes within 
the balls indicate their orientation. The parameters of the warp metric read R = 2. c = 1 , v = 2c. 



fy 



df(r) = ydf(r) 

dy r dr 



(A.8) 
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./; 



and 



dz r dr 

df(r) a[sech 2 (d(r + /?))-sech 2 ((j(r-^))] 



20 
(A.9) 

(A. 10) 



dr 2tanh(<7fl) 

Whether the velocity v of the warp bubble is time-dependent or not has only an influence on 
the Christoffel symbol Tf t , 



Appendix B. Geodesic equation 

The non-affinely parametrized geodesic equation reads ifTSIl 
c/V M dx v dxP _ 1 d£(o) dx» 
da 2 vp da da £(cr) da da ' 

where a is an arbitrary (non-affine) parameter and C, is a function of a. The affine parameter 
X then follows from 

X = / -yi—rda 
with 



da 

dX 






(B.l) 



(B.2) 



d 2 a 

dX 2 



1 dt, da 
a da dX ' 



(B.3) 



and an arbitrary but constant factor a. The relation between the affinely and the non-affinely 
parametrized geodesic equations can be derived by means of the chain rule 



dxP 

~dX 



dx^ da 
~da~dX' 



dV rfV (da 



dX 2 da 2 \dX 



dx» d 2 a 

da dX 2 ' 



(B.4) 



If we replace the non-affine parameter a by the coordinate time t, Eq. ( B. 1 1 with x° 
yields 

dx v dxP _ 1 d£(t) 



r° 

vp dt dt CO) dt 

Thus, the geodesic equation (|B.1[) can be written as 



d 2 x l 
~di 2 



-n 



dx v dx p 



-O 



dx v dx p dx' 



vp dt dt ' vp dt dt dt ' 



(B.5) 



(B.6) 



which simplifies to 



= 



d z x l 
dt 2 



dx' 
dt 



00" 



-2rt 



dxJ 



dx 1 dx k 



00" 



•2r; 



0j dt jk dt dt 

dx' ■ dx* dx k 

0j ~dT + jk ~d7~dT : 

where i — 1,2,3. The generalization of the constraint equation d9Jl, g^ v (dx^ /dX)(dx v /dX) 
Kc 2 with K = f or light-like and K = — 1 for time-like geodesies, reads 

dx' dx 1 



(B.7) 



dt \ ( dx' 

dX ( 800+280i ^ 



■Sij 



dt dt 



Kc~ 



The function C, follows from Eq. (B.5 1 



C(f) = Coexp 



i oo 



^o dx 
~ F °'d? 



o dx' dx' , 
ij dt' dt' 



(B.8) 



(B.9) 
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Appendix C. Euler-Lagrange equations 

The Euler-Lagrangian [16| equations for geodesies in the z = const hyperplane with 
Lagrangian 

^f = -c 2 i 2 + (x-vfi) 2 +y 2 (C.l) 

and v = const yield 

= Jx ^ + vf (i " V/ ^ " (i ~ vfi) ivft ' (C2) 

= ^j- [x - vfi] + (i - vfi) vif x , (C.3) 

0=y+(x-vfi)vif y , (C.4) 

where a dot means differentiation with respect to the affine parameter A . It is obvious that 
these equations are automatically fulfilled for 

i=h, x = vfk h y = k 2 , (C.5) 

with constants of motion k\ and k 2 . These initial values correspond to a null or time-like 
geodesic that starts perpendicular to the direction of motion with respect to the comoving 
reference frame of the center of the bubble. Because y grows linearly, the shape function / 
tends to zero and, thus, x is limited. Hence, such a geodesic approaches a line orthogonal to 
the x-axis. 

In case of a null geodesic, the initial direction reads k = ft) (±e(o) +C(2)), which yields 
k\ = ±co/c and k 2 = ft). Thus, there is no frequency shift. 

A time-like geodesic with initial local velocity Vp^ and four-velocity u = cy(e(o) + 

j8e(2)), where j3 = vp^fc and y = 1/yl — j3 2 , has constants of motion k\ = y and k 2 = y/3c. 
For null geodesies that are restricted to the x-axis, Eq. (C.4 1 simplifies to 

^f + v//=0 and j^i + v 2 // = 0, (C.6) 

where we made use of Jzf = 0, y = 0, and /, = — vf x . Thus, a future-directed initial direction 
k = ft) (em iem) leads toi — vi = ki = ft) [);(/ — 1 ) ± l] , which is a constant of motion, and 
we obtain 

i ± fc 3 . , C±vf 

t = 7Z -r, x = k 3 — . (C.7) 

cTv(l-f) cTv(l-f) 

Appendix D. Parallel transport of Sachs basis and Jacobi equation 

The parallel transport equation for the Sachs basis vector s = s^c^ can be cast into the form 



ds^ f n n dx' 

^i + { Qv+ iv Yt 

The Jacobian equation for the Jacobian field Y = Y^dn reads 

_ d 2 Y^ dY^ ( dxJ dx' dx 1 

" ~dli~ ~ ~dT ( r °° + 2r °' ~dJ + T ' k ~dT ~dt 



o = ~ + K v +r? v ~ )s v . (D.i) 



r M ,?rf dx' ^ dx' dx J \ v 
l 00,v+^0i,y dt +l U.v dt dt ) Y ■ 
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The initial Sachs basis vectors are perpendicular to the initial light direction k = —em) + 
cos^e^) +sin^e( 2 ), thus 

Si = — sin^ej!) +cos^e( 2 ) and S2 = e( 3 >. (D-3) 

1.2- 



The initial values for the two Jacobian fields read Y^ 2 \ . = and dY^/dt | „ = s 



Appendix E. Numerical integration 

The integration of the non-affinely parametrized geodesic equation avoids the problem of the 
inappropriate affine parameter. However, we still need the affine parameter for the calculation 
of the magnification factor. From Eqs. (B.2 1 and (B.5 1 we can expand the set of ordinary 
differential equations for X to 

dC, _( dx' dx l dx' ' 



-r°oo-24— -4 — — )C. (E.2) 



Since £ grows exponentially, which results in numerical problems, we substitute £ = exp(i//) 
in ( |E.2| i and obtain 

d\\f dx 1 dx 1 dx 1 

— = -r 00 - 2r 0( . — - rv. — — . (E .3) 

Here, we set a = 1 . 

Numerical problems due to exponential grow of the Sachs vectors and the Jacobi 
functions make the following substitutions necessary: 

s" = sinh (y ) , y = sinh («" ) . (E.4) 

The parallel transport of the Sachs vectors thus reads 

dp^ = ( ■ „ dx l \ mh{p°) 

<ft V ° CT to * / cosh(p" 

Note that there is no summation over the index jj. in the right hand side of this equation 



C+C— )— ^hk- (E.5) 



For the integration of the Jacobian equation (D.2 1, we obtain 
du^ 1 .. 



(E.6) 



— - = - tanh( M /J ){w ii ) (E.7) 

( u a dx'\ cosh(« v ) „ 

^u „^u dx* ^u dx' dx-*\ sinh(M v ) 



oo.v^ «,v A ' -y,v A A y cosh( M M)' 
As before, there is no summation over the index /i. 
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To apply numerical libraries like, for example, the Gnu Scientific Library [22] or the 
numerical recipes (23], we have to map the above equations to a one-dimensional array as 
follows: 

(E.8) 

(E.9) 

(E.10) 

(E.ll) 
(E.12) 

Here, n — 0, 1,2 and m = 0, 1,2,3. A light ray with initial direction k = — em) +cos^en) + 
sin^e(2) yields i = — 1/c, x = —vf/c + cosE,, and y = sin^. Thus y[3] — x/i = vf — ccosE,, 
y[4] = — csin^, and y[5] = 0. The corresponding Sachs vectors ( |D.3| l give y[6] = y[9] = 
0, y[7] = arsinh(-sin^), y[8] = arsinh(cos^), y[10] = y[ll] = y[12] = 0, and y[13] = 
arsinh(— 1). The initial values for the Jacobi vector fields are obvious. For the integration 
of the affine parameter A with A(0) = 0, we have £(0) = 1 and y/(0) = 0. Hence, y[30] = 
y[31]=0. 



r l 


= x" +l 


, y[» + 3] 




y[ n \ 


dt 


y[m + 6] 


— ri" 

- Pi ) 


y[m+10] = 


- n'" 

- Pi j 


y[ra+14] 


■vm 


y[m + 18] = 


dYl n 
~~dT' 


y[m + 22] 


- I \ ) 


y[m + 26] = 


dY{ n 
~~dT' 


ypo] 


= Y, 


y[31]=A. 





For the numerical integration of time-like geodesies, we can reduce the system (E.8 1 
to only eight equations. The array elements y[] can then be deduced from the following 
equations. Particles from the bridge have initial directions 

dt dx „ „ - dy „ . „ ,„ . „ N 

— =y, — =v/ + C j3cos^, -f= C /3sin^, (E.13) 

dk dt dt 

whereas particles injected from a static outside observer are described by 



dt 
In both cases, we have 1/(0) = In y. 



= yj3 \/c 2 — v 2 f 2 cos E, , — =c7j3sin^. (E.15) 

dt dt 
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